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ABSTRACT 

We use Suzaku observations to measure the spatial variation of the Fe Ka line with 
radius in the Tycho supernova remnant. The Fe line widths show a significant decrease 
from a FWHM value of 210 eV at the center to 130 eV at the rim. Over the same 
radial range the line center energy remains nearly constant. These observations are 
consistent with a scenario in which the shell of Fe-emitting ejecta in Tycho is expanding 
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at speeds of 2800-3350 km s . The minimum line width we measure is still a factor of 
two larger than expected from a single component plasma emission model. If thermal 
Doppler broadening is the dominant additional source of broadening, we infer an ion 
temperature of (1 3) x 10 10 K. 

Subject headings: ISM: individual (Tycho, SN 1572) — supernova remnants — X-rays: 
ISM 



1. Introduction 



Type la supernovae (SNe) play an important role in the chemical evolution of the universe by 
providing a significant fraction of the Fe group elements in stars, the interstellar medium and the 
intracluster medium. They have also become a prime tool to explore the expansion history of the 
universe. Yet, in spite of their importance, the physical processes involved in the actual explosions 
remain unclear. 

The SN seen by Tycho Brahe in 1572 can be classified with some degre e of confidence as Typ e 
la based on the light curve and color evolution from the historical record (jRuiz-Lapuente 1 12004 ) . 
The remnant (hereafter Tycho) is lo cated at Galactic coo r dinates (Lb) = (120° .0879, 1°.4460) and 



its estimated distance is 1.5-3.1 kpc (ICheyalier et al.lll 980: 



Ruiz-Lapuente II2004I ). but still debated (jSchwarz et al 



Albinson et al. 
1995h. 



1986 



Smith et al. 



1991 



Tycho is now a young and X-ray bright supernova remnant (SNR) that has been studied 
extensively for investigations of the SN la explosion mechanism. In the radio and X-ray bands, 
Tycho shows a relatively smooth and regular 8' diameter limb-brightened shell with several pro- 
tuberances and indentations notably toward the southeast (SE) and northeast (NE) parts of the 
shell, possibly driven by fingers of SN ejecta or the consequence of the blast wave interacting with 
a nonuniform ambient medium. The good angular resolution and photon statistics of Tycho ob- 
servations by Chan dra, and XMM-Newton have enabled detailed studies of the spatial structure 
(jWarren et al.ll2005l ) and have set significan t constraints on allo wed explosion models by compar- 
ing observed and simulated X-ray spectra (jBadenes et al.ll2006l ). Among the results revealed by 
Chandra is one that finds the region between the forward shock and contact discontinuity to be 
very narro w and dominated by nonthermal emission, concentrated in geometrically thin filamentar y 



structures (jHwang et al.ll2002l ; iBamba et al 



2005 



Warren et al 



2005 



Cassam-Chenai et al.ll2007l ). 



No thermal emission has yet been detected from the forward shock region; all the observed metal 
lines in the X-ray spectrum are produced predominantly by the ejecta. 



A direct measurement of the expansion velocity of Tycho has not yet been obtained. iReynoso et al 



(119971 ) reported the expansio n rate ob t ained from the radial displacement of the radio-emitting shell 
over the course of a decade. iHughed (|2000l ) obtained the rate from a difference of brightness pro- 
files measured from two ROSAT HRI observations taken in 1990 and 1995. These studies indicate 
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expansion velocities of the outer rim of Tycho of 2000-4300 km s _1 (radio) and 2200-4600 k m s - 1 
(X-ray), where the large range is dominated by the distance uncertainty. iRevnoso et al.l ([19971 ) 
found a variation with azimuth of the outer rim expansion rate and expansion parameter from 0.25 
to 0.75. Combining a distance of 2.3 kpc with their age (416 yr) and their radii leads to a range of 
observed rim velocities versus azimuth of 1700-4300 km s _1 . 

Tycho shows a strong Balmer line with a profile consisting of narrow and broad components. 
The width of the broad compon ent and intensity ratio have been used to measure the shock velocit y 
without distance dependence. ( Chevalier et al. 1980 ; Smith et al. 1991 ; Ghavamian et aDboOll ). 
These studies yield velocities of 1800-1900 km s _1 (H a width) and 1800-3500 km s _1 for the optical 
bright knot " g" . 

In this Letter we determine, for the first time, the expansion velocity of the Fe ejecta in 
Tycho from Doppler broadening of the Fe K lines, utilizing the modest imaging cap ability and well 



calibrate d instrument response o f the Suzaku CCD cameras (jOzawa et al. II2009I ). We refer the 



reader to lTamagawa et al.l (|2009i ) for more details on the broadband continuum emission of Tycho 
and the discovery of faint lines from the Fe-group elements Cr and Mn in the integrated Suzaku 
spectra. 



Observation 



Suzaku, the fifth Japanese X-ray astronomy satellite (jMitsuda et al.1 120071 ). was launched into 
a 550 km high orbit on 2005 July 10. It contains two types of functioning instruments: X-ray 
CCD cameras ( Kovama et al. 20071 , XIS) covering the 0.2-12 keV energy range, and a nonimaging 



hard X-ray detector (jTakahashi et al 



2007 



Kokubun et al.l 120071 . HXD) sensitive to X-rays in the 



10-70 keV band (PIN diodes) and 30-600 keV band (GSO scin tillators) . The four XIS cameras are 
each located at the focus of an independent X-ray telescope (jSerlemitsos et al.l 120071 . XRT). One 
of the four XIS cameras utilizes a back-illuminated (BI) CCD, whereas the other three use front- 
illuminated (FI) CCDs. (N.B., one of the FI CCD cameras mentioned was lost in 2006 November.) 

Tycho was observed by Suzaku on 2006 June 27-29 during the time allocated to the Science 
Working Group. The XIS was operate d in normal full-frame clocking mode without spaced-row 



charge injection (jUchivama et al. Il2009l ). 



We started with data processed through the standard pipeline (version 2.0.6.13) using calibra- 
tion reference files publicly released on 2008 February 1. Data accumulated during passage through, 
and for an additional 436 s after emerging from, the South Atlantic Anomaly (SAA) were discarded 
as were orbital periods of low cut-off rigidity (less than 6 GV). We also rejected data taken at low 
elevation angles to the Earth's rim (less than 20° from the limb of the sunlit Earth and less than 
5° from the rim of the dark Earth). Hot and flickering pixels were removed with SISCLEAN and 
only events with grades of 0, 2, 3, 4 and 6 were retained. 
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Suzaku also observed an offset region close to Tycho on 2006 June 29-30 for purposes of 
background estimation. We applied the same screening criteria and event selection described above 
to this background observation. Two faint unidentified X-ray sources appear in t he offset pointing 



rough ly 4' and 7' southwest of the center of the field of view (FOV) (see Figure 2 in lTamagawa et al 
2009). These sources are located in the Tycho annular analysis regions numbered 6 and 7 (see 



Section [3] below) . We can, however, safely neglect these sources, since their emission is 10 times 
fainter than that of Tycho in the 2-5 keV band. 

The net exposures of on-source and offset data after the data re ductions are 101 ks and 51 ks, 
respectively. Details of these observations have been summarized in lTamagawa et ajj (|2009l ). 



3. Spatially Resolved Line Spectra 



Guided by the nearly circular, shell-like structure of Tycho we divided the XIS FOV into a cen- 
tral region 1'.41 in radius, surrounded by six additional annular regions each 0'.54 wide in order to 
examine radial variations in the spectral parameters. Since the structure of the r emnant is distorted 



towa r d the SE, as c l early seen in the Chandra and XMM-JVewton observations (jDecourchelle et al 



2001 



Warren et al.ll2005l ). we have excluded the SE quadrant, as indicated in Figure [TJ The re- 



main ing NE, NW, and SW quadrants of Tycho show a similar shell-like morphology in the Fe K 
line (IWarren et al.ll2005l ). Background spectra were accumulated from the offset observation using 
the same region in detector space. 

Our spectral model consists of a power-law compo nent (which which well approximates the 
continuum emission in Tycho see lTamagawa et al.ll2009i ) and several Gaussians to account for the 
prominent emission lines. At the energy resolution of the XIS, most observed line features include 
multiple emission lines. In order to reduce the effects of line blending from different transitions 
or even different elemental species, the spectral fits were performed in the following narrow bands 
separately: 1.784-1.894 keV for He-like Si Ka, 2.40-2.60 keV for He-like S Ka and 6.00-6.90 keV 
for Fe Ka. For this Letter, we focus primarily on the Fe K line profiles. The lower energy lines, 
and Si He-like Ka in particular, are subject to complexities in the calibration of the energy scale, 
and thus require a much more careful analysis. We defer this to a subsequent paper, but use the 
Si and S line profiles for illustrative purposes only in this Letter. 

Best-fit values and uncertainties (at 90% confidence) are quoted in TableQ]for the Fe Ka spectra 
(labeled from the center, Regl, to the rim, Reg7) and Figure [2] shows the radial dependence of these 
parameters. The average centr oid energy we find, 6455 eV i s consistent with prey io us measurements 



by AS CA and XMM-Newton (jHwang and Gotthelfi Il997l : Badenes et al.ll2006l ). iTamagawa et al 



(J2009J) used the same Suzaku observation data and obtained the centroid over the whole remnant 
as 6445 eV. The difference of 10 eV is caused by using differe nt calibration reference files, since 
calibration reference files publicly released in 2007 were used in lTamagawa et al.1 (120091 ) . The line 
centroids prefer a slight radial gradient corresponding to an increase of ~12 eV from the center to 
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the rim. The line widths on the other hand show a highly significant gradient from a FWHM value 
of 210 eV near the center to values of ~140 eV at the rim. A flat profile can be rejected with high 
confidence (greater than 3a). 

To estimate the calibration error of the current response function, the width of the Mn Ka 
calibration line from the built-in calibration source ( 55 Fe), which illuminates the two far-end corners 
from the readout node, has been examined. The width of the calibration line has been obtained 
to be 16±i§ eV FWHM. In the source detected region, the broadening due to the charge transfer 
inefficiency (CTI) is smaller than that of the calibration source because of the reduced number of 
charge transfers. Our obtained widths (130-210 eV) in Tycho are significantly larger than that of 
CTI induced broadening. 

We note that because of the broad point-spread function (PSF) of the Suzaku XRT (half- 
power diameter of 2') the radial bins used here are not fully independent. To assess the effects of 
smearing, we used XISSIM to simulate an XIS image assuming a thin shell model f or the radial 
surface brightness profile of the Fe Ka line intensity (as given in lWarren et al.ll2005l ). Only 25% 



of the observed events in region 6 actually originate from this sky region; 60% come from regions 
4 and 5 where the Fe Ka brightness peaks. Virtually all of the photons in region 7 (the rim) come 
from regions 4, 5 and 6. Thus any true intrinsic spatial variation of the line parameters in the 
outermost three radial bins has been smeared out by the XRT PSF. On the other hand, smearing 
effects of the XRT are much less important for the variation interior to radii of ~3'. 

The ionization age can be estimated from the line center energy. Our measured central energies 
of 6448-6461 eV and widths of 130-210 eV FWHM indicate the existence of low ionized Fe below 
Fexvin- Of the nonequilibrium plasma emission models in XSPEC, only the NEI version 1.1 
model includes low ionized (below He-like) ions. We therefore use this model. Figure [3] shows the 
ionization age dependence of the centroid energy and width at various electron temperatures when 
the modeled Fe Ka blend is fitted by a single Gaussian line. 

Our measured central energies of 6448-6461 eV correspond to Fe ionization ages in the range 
10 10 < n e i/(cm -3 s) < 10 11 for all electron temperature values kT c > 0.6 keV. Simulated spectra 
over this range of thermodynamic parameters yield expected Gaussian widths for the Fe Ka blend 
of only 50-110 eV FWHM. These values are far below what we measure in the central region (~ 210 
eV) and even at the rim (>130 eV), implying the presence of additional line broadening mechanisms 
in TychoResults from the independent Si and S band analyses also show similar tendencies in 
the radial dependence of line centroids and widths. Calibration uncertainties and line blending 
complicate the analysis of these lower energy lines; we defer a quantitative discussion of these 
results to a future study. 

In order to directly compare the line profiles without model fitting, we took the ratio of the 
narrowband spectrum in region 7 (rim) to that in region 1 (center) for the Si, S and Fe Ka 
lines (Figure H]). All three ratios show a similar profile characterized by an "M" shape (lower in 
the middle, higher at the edges). The ratios show that the line shape in the central region is 
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significantly broader than at the rim for all three lines. 

To investigate further the broad Fe Kq line profile, we next fitted a double Gaussian line 
plus power-law continuum model to the spectrum extracted from the remnant's center. Both 
Gaussian widths were fixed at the value obtained from the rim (130 eV FWHM). The fits were 
acceptable as shown in Figure [5] with best-fit center energy values of 6405+^ eV and 6536±f| 
eV and normalizations of 3-llto 5 x 10~ 4 photons cm -2 s _1 and 2.1^q^ x 10 -4 photons cm -2 s _1 , 
respectively. The difference in line energies corresponds to a velocity difference of ±3040^240 km 
s _1 or it can be interpreted as contributions from ionization timescales of n e t ~ 1 x 10 9 cm~ 3 s 

and n e t ~ (5 10) x 10 10 cm -3 s. When we change the fixed width within the 90% confidence 

level (100-150 eV FWHM), the center energy values vary by only ±7 eV, well within the 90% error 
range. 



Discussion and Conclusion 



The excellent energy resolution of the Suzaku XIS has enabled us to discover significant broad- 
ening of the Si, S and Fe emission lines in Tycho as well as a clear decrease in the line widths from 
the center toward the rim. In the following discussion we consider, in turn, two simple astrophys- 
ical interpretations of these observations: (1) radial variations of the thermodynamic state of the 
X-ray-emitting plasma and (2) expansion of a spherical shell. 

The presence of two Gaussian lines in the central region can be explained naturally by the 
presence of two (or more) different plasma states: hot gas with low (n e t ~ 1 x 10 9 cm -3 s) and high 

(n e t ~ (5 10) x 10 10 cm~ 3 s) ionization timescales. There is some independent evidence for this, 

resting on the observation that the Fe Kg emission in Tycho peaks interior to that of Fe L and 



Si Ka (IHwang and G otthelf 



1997 



Decourchelle et al 



2001 



Hwang et al.l 12003 ). The theoretical 



interpretation (jBadenes et al.l 120061 ) is that the shocked ejecta exhibit a radial temperature and 
ionization timescale variation in the sense that near the contact discontinuity the ejecta are cool 
with high ionization age, while near the reverse shock they are hot with low ionization age. Guided 
by this theoretical picture, let us consider a simple two-zone model in which the interior half of 
the Fe-emitting shell produces a 6.41 keV line while the exterior half produces a line at 6.54 keV. 
The relative proportions of the emission from each of these shells seen in projection will vary from 
the center to the rim, which will produce a radial change in the observed line width. However our 
calculations show that the broad PSF of the telescope smears out the radial variation in width 
which, under this scenario, occurs most rapidly near the rim where the differences in the projected 
surface brightnesses of the emission zones are most extreme. This is unlike the smooth gradual 
drop with radius in the observed profile. Another problem for this model is that it produces a fairly 
large shift (greater than 30 eV) in the centroid of the Fe Ka line between the center and the rim. 
This is somewhat larger than the shallow gradient of ~12 eV observed in Tycholn Figure [21 profiles 
expected for the two-zone model (with assumptions of low ionized Fe shell (6.41 keV line) of 3'-3'.5 
and high ionized shell (6.54 keV line) of 3'. 5-4') are shown. It suggests that the two-zone model 
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cannot describe the constant center energy in the outer portion and gradual drop of the width. 

If the Fe emission comes from a radially expanding shell, then the projection of redshifted and 
blueshifted lines from the rear (receding) and front (approaching) parts of the shell will produce a 
double line near the center. Furthermore, the energy separation of the lines depends on the velocity 
component along the line of sight, which varies as a cosine function with projected distance from 
the center. At the rim the shell velocity is transverse to the line of sight and so the line width 
there should be minimum. In the ideal case the energy centroid of the line should be constant 
with radius, although deviations from spherical symmetry can introduce some variation. In general 
the centroid and width profiles of the Fe Ka line we measure in Tycho are more consistent with 
an expanding shell interpretation than multiple thermodynamic emission zones. In Figure [2j the 
radial profiles expected for the expanding shell model (with assumptions of expansion velocity of 
3000 km s -1 and radius of the shell of 3') are shown. 

The minimum line width of 130 eV FWHM we obtain at the rim, where velocity broadening 
should be minimal, is still wider than that expected from a single-component NEI model (50-110 eV 
FWHM). An additional source of broadening, at the level of 70-120 eV, is therefore required. Under 
the assumption that this comes fully from thermal Doppler broadening (i.e., ignoring macroscopic 
turbulence), we estimate the ion temperature to be kT lon = (1-3) xlO 10 K. This is consistent with 
the shock t emperature expecte d in models that well describe the ejecta spectrum of Tycho (see 



Our observations thus indicate that the Fe-emitting ejecta of the Tycho SNR is expanding 
at a velocity of ~3000 km s _1 . If macroscopic mixing of the ejecta has not taken place, then we 
would expect the intermediate mass elements (Si and S) to show even larger velocities. Estimating 
the expansion speed for these other species is going to be important for further investigation of 
the three-dimensional structure and dynamics of TychoFor that purpose, careful calibration of the 
energy scale around the Si Ka emission line and more precise modeling of the X-ray spectrum are 
needed. More detailed results including these other elements will be presented in a forthcoming 
paper (A. Hayato et al. 2009, in preparation). 

We thank all members of the Suzaku team for their careful work operating the satellite and 
calibrating the instruments. J.P.H. acknowledges support from NASA grant NNG05GP87G. A.H. 
and M.O. are financially supported by the Japan Society for the Promotion of Science. 
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Table 1. Best-fit single Gaussian fits to Tycho Fe Ka line 



Region Centroid Width (FWHM) 
(eV) (cV) 



Regl 


6456^5 


21 2+ 16 


Reg2 


6448±4 


186+.^ 


Reg3 


6452+J 


153112 


Reg4 


6454+J 




Reg5 


6457 + | 


132111 


Reg6 


6461+4 


146+| 


Reg7 


6458+4 


125_28 



Note. - 



All errors are quoted at the 90% confidence level. 



) 0.5 1 1.5 2 2.5 3 

Fig. 1. — 6-7 keV Fe-band X-ray image of the Tycho SNR obtained with the Suzaku XIS. Regions 
defined for spectral analysis (Section [3]) are indicated. 
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Fig. 2. — Radial dependence of best-fit line energies (a) and widths (b). Data points were obtained 
from spectra extracted from the sum of the NE, NW, and SW segments indicated in Figure [TJ Solid 
and dashed gray lines represent the roughly estimated profiles of the expanding shell model and 
the two-zone model, respectively. 
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Fig. 3. — Fitted Gaussian center energies and widths as a function of ionization timescale obtained 
from simulated nonequilibrium ionization plasma spectra with NEI model (version 1.1) at the 
electron temperatures of 0.6 (solid), 1.0 (dashed), 2.0 (dashed-dotted), 3.0 (dotted), and 4.0 keV 
(dashed-dotted-dotted-dotted). Horizontal lines in the upper panel show the range of observed 
center energies. The line in the bottom panel shows the measured FWHM of Region 5 for reference. 




Fig. 4. — Spectra and spectral ratios for bands containing the Si, S and Fe Ka lines. Top: spectra in 
regions 1 (open circles) and 7 (open squares). Spectra in region 7 were multiplied by 0.6. Bottom: 
ratios of spectra in region 1 to those in region 7. 
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Fig. 5. — Spectrum covering the 6.0-6.9 keV band with the best-fit double Gaussian and power-law 
model. Residuals (bottom panel) are in units of %. 



